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Using an iteration technique, we obtain exact expressions for the free energy and 
the magnetization of an Ising model on a two - layer Bethe lattice with intralayer 
coupling constants Ji and J2 for the first and the second layer, respectively, and 
interlayer coupling constant J3 between the two layers; the Ising spins also couple 
with external magnetic fields, which are different in the two layers. We obtain exact 
phase diagrams for the system and find that when \J^\ 0, ATc = '^"^'^^^jj^"^^^ ~ 
{Js/ Ji\^^^ , where Tc{J3) is the phase transition temperature for the system with 
interlayer coupling constant J3 and the shift exponent ip is 1 for Ji = J2 and is 0.5 
for Ji 7^ J2. Such results are consistent with predictions of a scaling theory. We also 
derive equations for ATc when IJ3I approaches 00. 



PACS numbers: 05.50+q, 68.35Rh, 64.60Cn, 75.10H 



1 



1. INTRODUCTION 



The physical properties of various magnetic layered structures and superlattices have 
been intensely studied both experimentally and theoretically for reasons ranging from fun- 
damental investigations of phase transitions to technical problems encountered in thin-film 
magnets . Experimentally, sub-monolayer and monolayer films of ferromagnetic materials 
offer challenging opportunities to fabricate materials with various novel magnetic proper- 
ties, such as giant magnetoresistance, surface magnetic anisotropy, enhanced surface mag- 
netic moment and surface magneto-elastic coupling. On theoretical grounds, surface mag- 
netism has been treated within several different frameworks: mean-field approximations , 
effective- field theories 0, spin-fiuctuation theory 0], renormalization-group methods [§, 
two-site cluster approximations and Monte-Carlo techniques 0. Though each method 
has it own advantages, there all have limitations in treating film systems. Numerical tech- 
niques such as Monte Carlo method can provide very accurate results for properties of finite 
systems; however, they are computation-intensive and can be carried out only for relatively 
small system sizes. 

Since exact solutions for realistic layer systems on regular lattices are generally un- 
available, one relies on approximation schemes to obtain a qualitative picture of the phase 
diagram. In our approach, we replace the original two-layer regular lattice by a two-layer 
Bethe lattice with the same coordination number q as the original lattice. Once this approx- 
imation is made, we can solve our model exactly as in the case of one-layer Bethe lattice 
p. It is now widely recognized (see e.g. |]9HT2|) that in many cases solutions of a spin 
model on Bethe or generalized Bethe lattices are qualitatively better approximations for the 
regular lattice than solutions obtained by conventional mean-field theories, because of the 



presence of correlations (albeit weak ones) in the former |T3|] and the lack of correlations 
in the latter. It has also been found that phase diagrams of an Ising model on a Husimi 



tree (a Bethe-like lattice) with two-spin and three-spin interactions |TO| , pn]] closely match 
exact phase diagrams of an Ising model on a two-dimensional Kagome lattice with two-spin 
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and three-spin interactions [0]. Of course, our approximation also has hmitations: since 



correlations are weak |T3[ , it predicts a transition temperature which is higher than that for 
a regular lattice, and it is usually not reliable for predicting critical exponents. On the other 
hand, the Monte-Carlo method will be highly reliable for predicting critical exponents. But 
we believe that our approach correctly gives the general shape of the phase diagram. 

It is well known that the Bethe and Bethe-like lattices cannot be embedded in a finite- 



dimensional Euclidean space without distortion in their bond angles and lengths . On the 
other hand, it has been pointed out by Mosseri and Sadoc |16| that these structures can be 
considered as regular lattices of fixed bond angles and lengths if they are embedded in a two- 
dimensional space of constant negative curvature (the hyperbolic or Lobachevsky plane H2 
ITtI). The surface of negative curvature are now being introduced to describe some complex 
structure with large cells, formed by inorganic and organic materials, which can be considered 
as crystals of surface and films. Among them are the cubic crystalline structures formed 



by amphiphilic molecules in the presence of water [jT8|,|T9| and a magnetically coupled three 



dimensional (Terephthalato) manganese(II) network ^U|,^. Furthermore, non- Euclidean 



hyperbolic symmetries have even been found in hexagonal and cubic-close-packed Euclidean 
crystals [^. The local structural similarities between that materials and negatively curved 



Bethe lattice suggest that the two-layer Bethe lattice considered in the present paper can 
be related to some physical systems. 

In this paper, we use an iteration technique to obtain exact expressions for the free 
energy and the magnetization of an Ising model on a two - layer Bethe lattice with intralayer 
coupling constants Ji and J2 for the first and the second layer, respectively, and interlayer 
coupling constant J3 between two layers; the Ising spins also couple with external magnetic 
fields, which are different in the two layers. Allowing for a difference in these two fields is 



important, because they act in an opposing manner on the zero-field boundaries [23|. We 



obtain exact phase diagrams for the system and find that when | J3I ^ 0, ATc = t (0) ^ 
\Jz/ Ji\^^'^ -I where Tc{Js) is the phase transition temperature for the system with interlayer 
coupling constant J3 and the shift exponent ip is 1 for Ji = J2 and is 0.5 for Ji 7^ J2. Such 



results are consistent with predictions of a scaling theory [p4H26[]. We also derive equations 
for ATc when | J3I approach 00. 

The theoretical works on thin-layer systems before the present paper were less general. 
The system of coupled two-dimensional Ising planes on regular lattices, e.g. the square 
lattice, is not exactly soluble; however, it has been investigated by a variety of approximate 
methods. Ballentine used high-temperature series expansions to study the model with 



Ji = J2 = J3 This work was later extended by Allan |^ to films up to five layers and 



by Capehart and Fisher |^ to films up to ten layers. The two - layer system where the inter- 



layer coupling constant differs from the intra-layer coupling constant was studied by Abe |24 
in the context of a scaling theory valid in the limit of a weak interlayer coupling. The more 
general case in which Ji 7^ J2 has also received some attention. The most complete treatment 
was that of Oitmaa and Enting who combined mean- field theory, scaling theory and 
high-temperature expansions in a study of the two-layer model, and calculated the variation 
of the critical temperature, the layer magnetizations and the interlayer correlation function 



with J3. Recently, Ferrenberg and Landau P2[ considered the same two-layer problem using 



Monte Carlo simulations and mean-field theory. The Bethe lattice version of a thin film was 
first studied in PB[, in which the phase diagrams of coupled bilayers with Ji = J2 and zero 



external fields were obtained. 

More recently, there have been many investigations of the two - layer Ising model due to 



stimulation from experiments for ultra - thin magnetic films [pq ,p^-^. The experimental 
study of such systems has made significant advances in recent years due to improvements in 
surface-force apparatus and microscopic techniques. Many theoretical studies were devoted 
to the model with two exchange parameters, i.e. Js for spins on free surfaces and J for all 
other spins. However, in our opinion the exchange interaction (J3) between the surface and 
the second layer has an important influence on the surface magnetic order. Therefore, in 
the present work, we study an Ising model on a two-layer Bethe lattice with three exchange 
parameters Ji, J2 and J3 and in the presence of magnetic fields, which are different in the two 
layers. We will calculate the free energy of that system and investigate its critical properties. 



In previous papers p4| , ^ , the relations between the free energies of a spin-1 Ising model 



on Bethe and Cayley trees and of a multi-site Ising model on Husimi lattices and generalized 
Cayley trees are obtained. In that approach one obtains the free energy from recursion rela- 
tions, and equations for physical quantities by differentiation of the free energy functional 
with respect to external fields. Recently, an elegant and original method of computing the 
bulk free energy for any model on infinite Bethe and Husimi lattices was presented by Gu- 



jrati In the present paper, we extend this method to study an Ising model on a bilayer 
Bethe lattice. In particular, we will calculate exact phase diagrams and the shift exponent 
ip for the system. 



In 1970 - 1971, Abe p4| and Suzuki [^] used a scaling theory to predict that the 



shift exponent ip for the two-layer planar lattice Ising model is equal to the susceptibility 
exponent (7 = 1.75) of the two-dimensional Ising model. In 1992, Angelini et al. predicted 
that = 1-5 [0; in 1993, Lipowski and Suzuki predicted that ip = 2 [Q. Very recently, 



Horiguchi and Tsushima used a /3 - function approach |^ to obtain ip = 1.73 ±0.04 for 
the two - layer square lattice Ising model with Ji = J2, which is very close to the theoretically 
predicted value ip = ■y = 1.75. They also found that the shift exponent for the system 
with Ji 7^ J2 is 0.5 and explained this value in terms of a scaling theory. In 1998 Lipowski 



47| used a transfer-matrix mean-field approximation to calculate the shift exponent ip. He 
found tp = 1.79 for the system with Ji = J2 and ip = 0.501 for the system with Ji 7^ J2. 

In this paper, we consider both Ji = J2 and Ji 7^ J2 for the Ising model on a two-layer 
Bethe lattice and obtain exact values of ip, which are consistent with the predictions of the 
scaling theory. 

The outline of this paper is as follows. In Sec. II, we present the bilayer Ising model 
and discuss different types of ground states. In Sec. HI, we derive the exact free energy, 
equations of state and order parameters for the Ising model on the bilayer Bethe lattice. In 
Sec. IV, we investigate temperature dependence of the order parameters and discuss the 
phase diagrams. In Sec. V, we calculate the critical temperature in weak and strong vertical 
coupling regimes and obtain exactly the shift exponent ip. In Sec. VI, we conclude with 



some general remarks concerning our results. 

II. TWO-LAYER MODEL AND ITS GROUND STATES 



In this section, we consider an ultra-thin film composed of two atomic layers Gi and G2 
such that the lattice structures of Ci and G2 are identical and each of them have sites 
and coordination number q; the corresponding lattice sites in Gi and G2 are labeled by i 
and i', respectively, where 1 < < N and they are nearest neighbors (nn) to each other. 
The Ising Hamiltonian on such two-layer lattice can be written as 

- I3H ^ Ji J2 ^i^j + cri'C^j' + Y Si<^i' + hi^^Si + h2'Yai', (1) 

<ij> <i'j'> <ii'> i V 

where (5 = (kBT)^^ with being the Boltzmann constant and T being the temperature, Si 
and CTj' take values ±1, Ji and J2 are, respectively, the coupling constants of the exchange 
interaction between the pair of nn spins in the first and the second layer, J3 is the coupling 
constant between a spin in the first layer and its nn in the second layer, and hi and h2 are 
magnetic fields acting on spins in the first and the second layer, respectively. 

This model has three order parameters, two of which correspond to the thermal average 
of total spins of the first and the second layers, respectively, 

N 1 ^ 

TT^i = J?12<^i>^ = -T^ 51 < ^i' > ■ (2) 

i=l i'=l 

These order parameters can be defined by variation of the partition function with respect to 
hi and /i2- The total magnetization density m and the density of staggered magnetization 
77 are defined by 

m^^{mi + m2), 77 = ^(^1-7712). (3) 

The third order parameter corresponds to the interlayer spin-spin correlation function be- 
tween nn spins of adjacent layers 

1 ^ 

P=1;t II (< ^i'^i' > -<Si>< >). (4) 

i=i'=l 
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Before studying the temperature dependences of the order parameters, let us investigate 
the ground states of the model at T = analytically. The ground-state energy in units of 
I Ji I and in the absence of magnetic fields may be described by the following Hamiltonian 



E = - E 



(5) 



<plaq> Ll-^ll Kll 

Here the summation goes over all plaquettes and each plaquette consists of four nearest- 
neighbor pairs of the two-layer system with one pair, < ij >, on Gi, one pair, < i'j' >, on 
G2, and two pairs, < ii' > and < jf >, connecting Gi and G2- 

By comparing the values of E for different spin configurations, wc obtain the ground- 
state phase diagrams shown in Figs. 1(a) and 1(b) for Ji > and Ji < 0, respectively. We 
find five types of ground states with following values of the order parameters (m, 77, p): 

(/) m = ±l, 77 = 0, p = 0, 

(//) m = 0, ri = ±l, p = 0, 

(III) m = 0, rj^O, p = 1, 

{IV) m = 0, ?7 = 0, p=-l, 

{V) m = ±l/2, ri = ±1/2, p = 0. 

The coordinates (J2/I Ja/g] Ji|) of the multiphase points are: 

Ai (0, 0), Bi (-1, 1) and Ci ^ (-1, -1), for Ji > (6) 

and 

A2 (0, 0), B2 (1, 1) and C2 (1, -1), for Ji < 0. (7) 

Phase (I) represents the usual ferromagnetic ordering mi — 1712 [r] — 0). 

Phase (II) represents ferromagnetic ordering in Gi and G2, but magnetizations in Gi 
and G2 arc antiparallel, i.e. mi = — m2 and m = (interlayer ordering is antiferromagnetic 
type) . It is worthwhile to note that this phase corresponds to the well known compensation 



phenomenon which occurs when the magnetizations of two layers cancel each other instead 
of being equal. 

Phase (III) represents the antiferromagnetic ordering in both layers (mi = 1712 = 0) 
where interlayer ordering is ferromagnetic (p = 1). 

Phase (IV) represents the totally antiferromagnetic ordering (p = —1). 

Phase (V) represents the ferromagnetic ordering (m = ±1/2,1] = ±1/2), which is equiv- 
alent to the case that the ground state of one layer is ferromagnetic and the ground state of 
another layer is antiferromagnetic. 

The phases (I) - (V) will be referred to as (F) - ferromagnetic, (C) - compensated, (M) 
- mixed, (A) - antiferromagnetic and (SF) - surface ferromagnetic phase, respectively. 

III. EQUATIONS OF STATE AND FREE ENERGY. 

Let us consider an Ising model on a bilayer Bethe lattice, which is constructed by connect- 
ing to the central pair of sites q pairs in order to form the first generation and by connecting 
successively (g — 1) pairs to each pair in a generation to form the next generation. The 
result is an infinite lattice in which every site has (q + l) nearest neighbors, where q nearest 
neighbors are in the same layer as the site and one nearest neighbor is in another layer. 

The partition function of the system represented by Eq. (|l|) may be written as 

{(7,s} [ <ij> <i'j'> <ii'> i i' ) 

where the sum goes over all configurations of the system. 

Now we derive exact recursion relations for Z. When the Bethe tree is cut apart at the 
central pair, it separates into q identical branches, each of which contains {q — 1) branches. 
The partition function can be written as follows: 

Z = exp{J3aoSo + hiSo + h2ao}gl{ao,So), (9) 

{(TO, so} 

where ctq and sq are the spins of the central pair, n is the number of generations {n — >■ 
oo corresponds to the thermodynamic limit where surface effects may be neglected) and 
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(7n(o"0)So) is the partition function of a separate branch. Each branch, in turn, can be cut 
apart at the pair of sites nearest to the central pair. The expression for 5'n(o"o,So) can 
therefore be written in the following form 

gnicro,So)= exp{JiSoSi + J2croai + JsaiSi + hiSi + h2ai}g'J,Z\icri,Si). (10) 

{ci.si} 

Let us introduce the following variables x„, i/n and tn 

_ 9n{ + + ) _ gni + -) 

9n{ ) ^n(- + ) 

^From Eq.(|TO|) we easily obtain the recursion relations: 

■^n f\^^^n—\iyn—\itn—\)i 
Vn = /2(X„-1, Vn-l, tn-l), (H) 
tn = fsi •^n—li Vn—li tn—1 ) ; 

An exp ( Ji + J2) + exp (- Ji - J2) + -D„ exp ( Ji - J2) + -B„ exp (- Ji + J2) 
An exp (- Ji - J2) + exp ( Ji + J2) + exp (- Ji + J2) + -B„ exp ( Ji - J2) ' 

An exp (- Ji + J2) + exp ( Ji - J2) + -D„ exp (- Ji - J2) + -B„ exp ( Ji + J2) 
A„ exp ( Ji - J2) + exp (- Ji + J2) + -D„ exp ( Ji + J2) + i?„ exp (- Ji - J2) ' 

An exp ( Ji - J2) + exp (- Ji + J2) + -D„ exp ( Ji + J2) + -Bn exp (- Ji - J2) 
y4„ exp (- Ji - J2) + exp ( Ji + J2) + -D„ exp (- Ji + J2) + -Bn exp ( Ji - J2) ' 

An = xt' exp (2/ii + 2/12), = C exp (-2J3 + 2/ii), 5„, = ylHt^ exp (-2 J3 + 2/12). 

Through and tn one can express the magnetization and other thermodynamic 

quantities, so we can say that in the thermodynamic limit {n —>■ 00) Xn, Vn and tn determine 
the states of the system. For this reason the recursion relations can also be called the 



9n{- + ) 
9n{ )■ 



where 

fl(yXn,ynjtn) 

f2{Xn,yn,tn) = 
fsiXn^Un^tn) 

with 
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equations of state (EOS) for the two-layer Ising model. The magnetizations of the first and 
the second layers as well as the spin-spin correlation function between spins of adjacent 
layers are expressed by: 

^1 =< So >= =; —, (12) 

1712 =< O"0 >= Z 7^ ^1 (l>j) 

< CToSo >= -j =; —. (14) 

We are interested in the case when {xn,yn,tn) converges to a stable point {xs,ys,ts), 
which is associated with the thermodynamic solutions of the two - layer Ising model. In this 
case the recursion relations (or equations of state) given by Eq. (|11|) can be rewritten in the 
following form: 

fi^)'"exp(2/,,-2ft0 = ^. (15) 

\l+yj ui + vi 

i^)'-exp(2/..2/.,^^, (16) 
^^"-'(l^)'"e-^P(-«3)^|^. (17) 



where 



Ml = Cix — C2, Vi = sixy + S2t, Ci = cosh (Ji + J2), Si = sinh (Ji + J2), 
U2 = ci - C2X, V2 = Sit + S2xy, C2 = cosh ( Ji - J2), S2 = sinh ( Ji - J2), 



and 



— -I I ''SI y 1 I 1 



l + Xs 1 + 2/s l + Xs 
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The total magnetization density (m = (mi + m2)/2), the density of the staggered mag- 
netization [rj = (mi — 1712)12) and the density of the interlayer spin-spin correlation function 
(p =< (ToSo > — mim2) can be expressed as 

tU2 + f 2 



m = , (18) 

(19) 



U2 + tV2 + XUi + vixy'' 
[xyui + xvi) 



U2 + tv2 + xui + t;ia;?/ 

^ (l-t^)(»i-x;2^)-x^(l-t/^)(t.f-t;i^) _ .20) 
(u2 + tf2 + a;ui + xyf i)2 

In the case when (x„, t„) converge to a stable point (x^, y^, t^), we can obtain an equation 

for the free energy functional F: 



-f^F = -\ In {u\ - vDiul - vl) + ^ lna:^(l - t^)(l - y') 

q — 2 1 q 
— In {u2 + tv2 + xui + xyvi) + - ln2 + - In |ci - cl\. (21) 

In deriving this equation we have used the exact relation between the free energy of the 
Bethe lattice and Cayley trees |T2|,H3,ES . 



It is easily seen that the expressions for the order parameters mi, m2 and p can be 
obtained by differentiation of the free energy functional of Eq. ( ^I]) with respect to the 
magnetic fields hi, /i2 and the coupling constant J3, respectively. In this sense, the interlayer 
coupling constant, J3, is analogous to an external field. 

This result for the free energy is very useful for locating phase transitions in case of 
multiple solutions of the equation of state and for determining the equilibrium state. Using 
this free energy functional one can obtain the full phase diagram, describing not only the 
continuous phase transitions but also the discontinuous ones. 

In the next section we will discuss the critical properties of our model; in particular, we 
will calculate the critical temperature as a function of ratios of coupling constants and will 
show the full phase diagram in the three-dimensional parameter space spanned by coupling 
constants Ji, J2 and J3, for different values of the coordination number q. 
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IV. PHASE DIAGRAMS 



Now we consider the critical properties of the Ising model on a two-layer Bethe lattice 
with different ferromagnetic coupling constants (Ji > 0, J2 > 0). Without loss of generality, 
we need only consider Ji> J2- The phase transition occurs when hi = h2 = 0. In this case 
Eqs. (|15D, (|16D and become: 



X 



2(9-1) 



2\ 



exp (-4J3 



2 2 ' 



(22) 



(i + yy-\m-vi) = {i-yy-\m + vi) 



(23) 



(1 + ty-\u2 - V2) = (1 - ty-\u2 + V2) 



V2 = tU2^it'^), 



(24) 



where 



[^1 



[^1 



Hx) = I E Cl-lix- \/\Y: CVx^^ , , 

n=0 



n=0 



and C| ^ is a binomial coefficient. 



It can be seen that there always exists one solution of this system: 



y = t = 0, and x'^ ^ exp (— 2 J3 



Ui 
U2 



(25) 



This solution corresponds to the high temperature paramagnetic phase (mi = 0, m2 = 0). 
In addition to the first solution, the equations of state also have other solutions with y 7^ 
0, t^O: 

Ml$(?/2) - Six] \u2<^{t'^) - Si] , 



X 



2(q-l) 



U2<^{t' 



l-t2 



Si 



2\ 9-1 



xy 



ui^{y^) - six 



exp (-4J3) 



Of course, only the solution which minimizes the free energy functional (^Tj) is thermody- 
namically stable. The others correspond to unstable or metastable states. If there are two 
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or more solutions which have the same minimum free energies, these phases coexist and the 
system has a first-order phase transition. 

When these two solutions merge into one solution, i.e. {i) = {ii), we obtain the critical 
line of second-order phase transitions 

exp (2 J3) = ^^^^xr \ (26) 

CiXx - C2 

where xx is the solution of the following equation 

slxx = [((g - l)ci - si)xx - (g - l)c2][(g - l)ci - si - (g - 1)c2Xa]. (27) 
It is convenient now to introduce the new parameters ki and ^2 



ki = \ tanh Ji tanh 



2, 



k2 = Jtanh(-Ji + 0.51n^--)tanh(-J2 + 0.51n^--). (28) 
V Q ^ Q ^ 

The two solutions x^^''^^ of the Eq. (pTl) can thus be expressed as 

(1) 1 - kik2 (2) 1 + kik2 , . 

x\ = r-;— , x\— r-;— , (29) 

and the corresponding expressions for the A-lines of the second-order phase transition in the 
three-dimensional parameter space spanned by Ji, J2 and J3 will take the form 



q-l 



The critical lines of the second-order phase transition given by Eqs. ( |5(]| ) and (^Tj) separate 
the paramagnetic (P) phase from the ferromagnetic (F) and compensated (C) phases, which 
are in turn separated by a first-order phase transition line. 

Before discussing the phase diagram, it is convenient to introduce the parameters n = 
J2/J1, S = -h/lJi and T = Jf^. In terms of the T, n and 5, Eqs. ( pOD and (^TJ) for the A - 
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lines imply a relation T = Tc{n, S) which locates the critical temperature as a function of n 
and 6 for arbitrary values of the coordination number q. The two critical lines start at 

rpmax 2(1 -\- n) . 

=ln[g/(g-2)]' ^ ^' ^^'^ 

and meet each other at 



At J3 = the system has a second critical point 

„„„„ 2n 



h = 0. (34) 



^ ln[g/(g-2)]- 

The phase diagrams (T^ versus 5) of the Ising model on the two-layer Bethe lattice for 
different values n = 1,0.75,0.5,0.25,0.1 and for different values of coordination number 
q = 3,4,6,00 are shown in Figs. 2(a), 2(b), and 2(c). A few comments are in order. For 
J3 = we recover two critical temperatures of two single-layer Bethe lattices with different 
intralayer ferromagnetic coupling constants (Ji and J2). In the opposite limit of | J3I — > 00, 
the critical temperature goes asymptotically to a value given by Eq. (|52D with the effective 



intralayer coupling constant Ji(l + n), since the interlayer pairs become rigidly correlated. 
V. WEAK AND STRONG INTERLAYER COUPLING REGIMES 

The spin-1/2 Ising model on a two-layer square lattice is exactly soluble only in the cases 
J3 = and I J3I 00, where it is related to the one layer square Ising model. When J3 = 0, 
the Hamiltonian given by Eq. (0) describes two uncoupled Ising model or, equivalently, 
two free fermionic fields. In strong vertical interaction limits | J3I — > 00, each pair of spins 
coupled across the layers will act as a single spin, and the Hamiltonian given by Eq. (|I]) 
describes an one-layer Ising model with (Ji + J2) as the coupling constant. 

In the weak interlayer coupling regime (J3 ^ 0) the shift exponent ip can be defined by 

- T,(J3) -Te(0) /o.x 

Ar. = |J3Mh, (35) 
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where Tc{J^) is the critical temperature when the system has interlayer couphng constant 

In this section we calculate exactly the shift exponent for the Ising model on a two-layer 
Bethe lattice. In the weak coupling regime we obtain 

l^3| 



Ar, = 6i(g). 



and 



Ji 



for Ji = J2 {n = 1) 



(36) 



ATe = 62(g,n)(-^ 



for Ji > J2 (n < 1) 



(37) 



where 



hi{q) 



and b2{q,n) 



In a ( 1 + d 



n-l ' 



- (a^" - 1) 



g-2 "^^"^'-^^ 8(g- 1) \1 - a'- 

with a = g/(g — 2). 

Thus we find that the shift exponent ip for the system with Ji = J2 is equal to 1, which 
coincides with theoretically predicted results ■?/' = 7 = 1 for the single-layer Bethe lattice. 
For the system with Ji 7^ J2, we find that ip = 0.5, which also exactly coincides with the 
value predicted by the scaling theory ||2^ . 



In the strong coupling regime we have: 

TciJs) . ^ / 2IJ3MI 
1 — A exp 



c 



c 



for Ji > J2 {n< 1) 



(38) 



where 



with 



K 



2(g-i: 



ln[g/(g-2)] 



2{q-l) 



n + l 



n + 1 



+ 



It is easy to see from Eqs.(^), ( pTf ) and (^) that the behavior of the strong coupling 
expansions is very different from the behavior in the weak coupling regime. It seems that 
we are the first group to obtain Eq. (pSf ) for the two-layer system with different intralayer 
coupling constants (Ji 7^ J2). It should be noted that for the case Ji = J2, equations similar 
to Eq. (^) had been obtained by approximate methods fS^^^ . 
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VI. SUMMARY AND DISCUSSION 



In the present paper we have investigated an Ising model on a bilayer Bethe lattice with 
intralayer coupling constants Ji and J2 for the first and the second layers, respectively, and 
interlayer coupling constant J3 between the two layers. We first analyze phase diagrams 
of ground states, then using an iteration technique to obtain exact expressions for order 
parameters and the free energy of the bilayer Ising model (Eqs. (18)-(21)). We then obtain 
exact phase diagrams of Eqs. (30) and (31) and analyze these equations in the weak and 
strong interlayer coupling regimes, see Eqs. (36)- (38). The shift exponents in Eqs. (36) 
and (37) are the first exact result to support the scaling theory for ip, which states that ip is 
equal to the exponent of magnetic susceptibility for Ji = J2 and is equal to 0.5 for Ji 7^ J2 
p^^6[. It seems that Eq. (38) is a new result. 

In Sec. II, we present very rich phase diagrams for ground states. However, in Sec. IV 
we consider only phase diagrams for Ji > J2 > 0. It is of interest to study the evolution 
of phase diagrams in Sec. II as the temperature increases from to high temperatures. 
However, the analysis of such general phase diagrams is quite complicated. 

The dependence of various quantities on the film thickness is a topic of current interest. 
In principle, we can extend our calculations from two layers to n layers. For such a n-layer 
system, we can introduce n external magnetic fields hi, /i2, . . ., {hi for the i—th layer 
with 1 < i < n), 6*2 (= n{n — l)/2!) interlayer coupling constants for two-layer coupling, C3 
interlayer coupling constants for three-layer coupling, . . ., and C"(= 1) coupling constant 
for n-layer coupling. Therefore, the total number of such field coupling parameters are 

+ C2" + + . . . + Q = 2" - 1. Equations (15)-(17) for two-layer systems for three field 
coupling parameters {hi, /i2, ^3) can be extended to 2" — 1 equations for 2" — 1 field coupling 
parameters. It is very difficult to find analytic or numerical solutions of these equations for 
n > 2. However, we can simplify the problem by reducing the number of independent field 
coupling parameters, i.e. setting hi = h for 1 < i < n and keeping only a two-layer interlayer 
coupling parameter for two nearest-neighbor layers. We are working in this direction. 
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FIGURES 

FIG. 1. The ground-state phase diagram of the two-layer Ising model for (a) Ji > and (b) 
Ji < 0. 

FIG. 2. Phase diagram on {J^/qJi,T/q) plane for an Ising model on a two - layer Bethe lattice 
with intralayer coupling constants Ji and J2 for the first and the second layer, respectively, and 
interlayer coupling constant J3 between two layers; here q is the coordination number for one-layer 
Bethe lattice and T = 1/Ji > 0. (a) Ji = J2 and g = 3. A first-order phase boundary (dashed 
line) separates two ordered phases designated by (F) and (C). The solid line denotes the second 
order phase transition line, which separates paramagnetic phase (P) from two ordered phases (F) 
and (C). Notice that Tc is the critical temperature of the Ising model on one-layer Bethe lattice, 
(b) J\ = J2 and g=3, 4, 6, and 00. (c) g = 3 and n = J2/Ji=l, 0.75, 0.5, 0.25, and 0.1, which are 
denoted on curves by 1, 2, 3, 4, and 5, respectively. T^^ T4, and T5 are second critical point of 
Eq.(34) for curves labeled by 2, 3, 4, and 5, respectively. 
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